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ABSTRACT 

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galax- 
ies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright 
(—15 < My < —8) satellite galaxies. These simulations resolve high density regions, comparable to 
giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H2 
formation and destruction that ties star formation to the presence of shielded, molecular gas. Before 
infall, supernova feedback from the clumpy, bursty star formation captured by this physically moti- 
vated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the 
massive satellite progenitors {M^r > IO^Mq, M* > IQ'^Mq) compared to DM-only simulations. The 
progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to 
both heating at reionization and gas loss from initial star forming events, preserving the steep inner 
density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce 
the central densities of the luminous satellites, particularly those that enter with cored dark matter 
halos, increasing the discrepancy in the central masses predicted by baryon-f DM and DM-only simu- 
lations. We show that DM-only simulations, which neglect the baryonic effects described in this work, 
produce denser satellites with larger central velocities. We provide a simple correction to the central 
DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations 
should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf 
satellites. 

Subject headings: Galaxy — halo; galaxies — dwarfs 



1. INTRODUCTION 



The favored Cold Dark Matter (CDM) cosmologi- 
cal model has been successful in reproducing many 
large scale observable properties of the Universe (e.g., 
Efstathiou et al.l 119921 : IRiess et al.l 119981 : ISpergel et atl 



20071) . The CDM model, however, still faces many 



challenges from observations of galaxies on small scales. 
The most well known of these problems has been 
termed the "missing satellite problem", since CDM- 
based models predict orders of magnitude more dark 
matter subhalos within the virial radii of Milky Way- 
massed galaxies than are observed as luminous satel- 
lites of such systems (iMoqre et afl 119991 : iKlvpin et al.l 
119991: IWadepuiil fc Springell 120111) . Another aspect of 
the missing satellites problem is the discrepancy between 
the masses of the most massive pred icted subhalos and 
the most ma s sive observed satellites (iMoqre et alJll999l : 
IKlvpin et al.l [1991 iBovlan-Kolchin et al.l I2011L I2012D . 
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Simulated ^ IO^^Mq halos consistently have several sub- 
halos that are too massive and too dense to host the most 
lumin ous dwarf spheroidal (dSph) satellites of the Milky 
Way fBovlan-Kolchin et al.! 12011 1201I IWolf fc BullocH 



[2012: Hayashi & Chiba 20121). This appears to provide a 
new challenge for the CDM paradigm on small scales, be- 
cause the most massive subhalos of an L* galaxy should 
not be devoid of stars. 

The tension between the predicted and observed inner 
densities of the Milky Way's dSph satellites is reminiscent 
of the longstanding tension between the predicted and 
observed shapes of the central density profiles of galaxies 
(known as the "core/cusp problem"). The steep inner 
density profiles and concentration s predicted for dark 
matter halos and their satellites (iNavarro et al.l 1997 ; 



Stoehr et al.|[2q0 ^iDck el et al.ll2003bll a': 'Reed et al.lbooa 
Springel et 311120 08: Madau et al.l 12008: Na varro et al.l 
2OIOI : iMaccio et al.l |2009|) are inconsistent with those 



obser v ed in isolated field ga la xies (Ivan den Bosch et al.l 
2OOOI: Ide Blok et all [2OOII: Ide Blok fc Bosmal 120021: 



Simon et al.l 20031 iS water s et al. 1 120031: IWeldrake et al 



2003 : 
2007 



iKuzio de Na rav et al 
lSpanqetal.1 120081: 



[200a 



deBloket al.1 120081 : lOh et al.l 
2OI2D . and observed in s atellite s (iKlevna et al 



2011b IDel 



Gentile et al. 
Trachternach et al.l 12001 



iMashchenko et al.l [20051 iGoerdt et al 



opolql 



Strigari et all 120061: 



20091: iStrigari et ah! 



2002 

2ooe 



Jardel fc Gebhardtl I2012t " 
Havashi fc Chibal I2012t) 



Gilmore et al.l l2007t iWalker et al.l 
2010t iWalker fc Pefjarrubial 



2011 



2012 



IWolf fc Bullockl 
This inconsistency is in- 
dependent of whether the density slope in simulations 
follows log(p) oc 7 log(radius), with 1.0 < 7 < 1.5 (i.e., a 
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NFW profile, iNavarro et allll997[). or a pow er-law slope 
(i.e., an Einasto profile. iNavarro et al.ll20l"ol ). 

Baryonic processes have often been proposed to ad- 
dress the apparent discrepancies between observations 
and the predictions of DM-only simulations. Within 
the CDM paradigm, the missing satellite problem is 
likely reconciled through reionization, which is expected 
to suppress star formation in subhalos wi th maximum 
circu la r velocity, Vr,, at in f all < 30 km/s (iQuinn et al 



1996; 



1997 



iThoul fc Weinberg' '1996" 
Gnedin 2000: Hoeft et al 
20081: IMadau et al.Tl200ah . and supernova feedback fur 



''Navarro & Steinmctz 
2006: Okamoto ct al 



ther suppressi ng star form a tion at the mo re massive 
subhalo end (iDekel fc Silkl [l98l iB cnson c t~all [20021: 
Dekel fc Wo"oll2003t iGovernato et aL[i2007i: iBusha et al.l 



201 nfi . When combined with observational incomplete- 



^ Q Piete- 

ness 'effect s ('Simo n fc Gehal [2007[ iTollerud et al.l [2001 
IWalsh et al. 2009: .Koposov et al.ll2009( ). it is likelv that 



baryonic effects can bring the predicted number of sub 
halos in line with the observed number of satellites. Yet 
even if the correct number of satellites can be repro- 
duced, CDM still predicts that the most massive satel- 
lites today are more massive and more dense than ob- 
served for Milky Way satellites (iBovlaii -Kolchi n et al.l 
[20T1 iWolf fc Bullockll20ll iHavashi fc Ctiiba.20121) ! 

Baryons have also been invoked to reconcile the 
cusp/core problem in CDM. The predicted cuspy pro- 
files of dark matter halos may not be physically realistic 
because the effect of energetic feedb ack from supe r novae 
(SNe) can reduce t h e baryonic (iDekel fc Silkl 119861: 
Mailer fc Dekel [200l IDekel et al.l l2003al: iBrook et al.l 



20111: iGuedes et al.l I2011D and DM mass at the cen 



ters o f~Kalaxies (INavarro et al.l 1996 : Read fc Gilmore 



20051: iMashchcnko et al 



20101 : ^Pasctto e t al.l 



Cloet-Osselaer et al.l 



Pontzen fc Governatol 




2008: Governato et al 



de Souza et al.l 2011 



iMaccio et al.l 



iGoyernato et al.l [2012 



2012 



Tevssier et al.ll2012l: lOgiva fc Morill2012D . Importantlv. 



there is not yet a unifying baryonic solution that solves 
the missing satellites, cusp/core, and massive subhalo 
problems simultaneously. 

In fact, unlike the missing satellites problem and the 
core/cusp problem, solutions to reconcile the predicted 
massive subhalos with observed satellites have not yet 
addressed the effects of baryons. For example, studies 
have concluded that the excess of predicted massive sub- 
halos around Milky Way-mass galaxies may disappear if 
satellites are modelled with Einasto (versus NFW) den- 
sity profiles, or if the Milky Way's tru e virial mass is 
8 X 10" Mc, rather th an 10^^ Mq (|Vera-Ciro et al.l 



120121: [Wang et al] 120121 ). It is not yet clear whether 
these solutions are themselves sufficient to fully reconcile 
observations with CDM-based models, and alternative 
cosmological models (e.g., warm or self- interacting dark 
matter) have also b e en proposed to explain observations 
(IAlam etal.1 [20021: iStrigari et al.l l2007bl: iLovell et al " 



2012: iVogelsberger et al.l [20121: Ivan den Aarssen et~ 



20121 ) 



The inclusion of baryonic physics thus remains a glar- 
ing gap in predictions for the properties of dwarf satel- 
lites. Here we investigate whether the inclusion of en- 
ergetic feedback from stars and supernovae, which has 
been shown to reduce the central dark matter densi- 
ties of simulated field galaxies, can also reduce the pre- 



dicted central densities of luminous satellites around 
Milky Way-mass galaxies. One reason that baryonic ef- 
fects have not been studied in detail in satellite galaxies 
is because including baryons in cosmological simulations 
makes them much more c omputatio nally expensive. Pre- 
vious work (IGovernato et al. 2010t iGuedes et al.l 120111 : 
[Pontzen fc Governatol l2012h demonstrated that simply 
achieving high resolution is not enough to reduce the cen- 
tral concentration of dark matter in galaxies. Instead, it 
is necessary to also adopt a physically motivated model 
in which star formation is limited to high density peaks, 
with densities comparable to the average densities in a 
giant molecular cloud. This model allows highly over- 
pressurized regions to form when supernova energy is in- 
troduced to the high density surrounding gas, and the 
resulting hot bubble of gas fiattcns the central potential, 
leading to irreversible e xpansion of dark matter orb its as 
well as driving a wind ([Pontzen fc Governatol[20T^ . 

The implications of such feedback may be all the more 
important when predicting properties of the satellites of 
larger galaxies (such as the Milky Way's own dwarf satel- 
lite population) , because any effect imprinted on the cen- 
tral densities of the dwarf galaxies by baryons could be 
exacerbated during the dwarf galaxies' tidal evolution 
around their eventual host. A proper prediction for both 
the total number and the internal structure of the Milky 
Way and MSl's dwarf galaxies may thus require an ac- 
curate treatment of baryonic feedback. Only a few DM 
+ baryon studies have focused on the satellite popula- 
tion of a Milky Way-massed galaxy, but at much lower 
resolutions than have been achieved in DM-only simula- 
tions, making comparisons with DM-only predictions for 
inner density pr ofiles particularly difficult (see, however, 
iDi Cintio etal][20Tl(Parrv et al.llMl ). 

In this paper, we compare DM-only and Smooth Par- 
ticle Hydrodynamics (SPH) simulations of two Milky 
Way-massed galaxies run in a cosmological context to 
z = 0. In fact, the simulations used in this work 
adopt a model that has already been shown to al- 
leviate some of the small scale problems of CDM 
(IGoyernato et al.l 120101: lOh et all boiH IBrook et al 



20111: IMaccio et al.l 120121 : [Pontzen fc GovernatoTIm-, 
Governato et al.ll2012D . This model is capable of forming 



i 



bulgeless dwarf disk galaxies that have "cored" dark mat- 
ter density profiles (i.e., shallower inner density slopes 
than predicted by CD M simulations without baryons, 
IGovernato et al.ll2010L [2012: Oh et al. 20lJ[). While the 
previous studies have focused on isolated field galaxies, 
this paper examines how this same model affects the 
satellite population of Milky Way-mass galaxies. 

This paper is organized as follows. In Section 2 we de- 
scribe the simulations used in this work, as well as our 
method of selecting the satellite sample. In Section 3 we 
focus on the z — density and circular velocity profiles of 
SPH satellites, and compare to their DM-only counter- 
parts. Section 4 explores the role of supernovae feedback 
and tidal stripping on the mass dependent evolution of 
satellites, from high redshift to z = 0. In Section 5 we 
propose an update to the standard treatment of satel- 
lites in DM-only models. We summarize our results and 
conclude in Section 6. 



2. SIMULATIONS 
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The high-resolution N-body + SPH si mulations used 
in th is paper were run with gasoline (jWadslev et alj 
l2004f ). The two halos studied in this paper were ini- 
tially selected from a uniform resolution, DM-only, 50 
comoving Mpc box. The initial conditions for this box 
were generated ass uming a WMAP Year 3 cosmology 
(|Spergel et alJ[2007l) : ilra = 0.24, = 0.76, Hq = 73 km 
s-i, erg = 0.77. The two halos (h277, M^.^^ = 7xlO"Mo, 
and h258, M^ir 8 x IO^Mq) were selected based on 
their z = mass and merging histories. We define virial 
mass relative to cr itical density , pc, using p/pc = 100 
at 2 = following [Grosd (jl997T) . Galaxy h277 was se- 
lected to have a quiescent merger history, with its last 
major merger occurring at z ^^3, while h258 has a binary 
merger aX z = \. Each of h277 and h258 was resimulatcd 
at higher resolution (and with ga s particles) using th e 
volume renormalization technique (jKatz fc Whitelll993D . 
This approach simulates only the region within a few 
virial radii of the primary halo at the highest resolution, 
while still maintaining the large volume at low resolu- 
tion in order to account for the large scale tidal field 
that builds angular momentum in tidal torque theory 
(|Peeblesll 19691 iBarnes fc Efstathioul[T987h . The high res- 
olution simulations were run from = 150 to z = 0. 

N-body-|-SPH volume renormalized simulations of 
both h277 and h258 have been studied in previously pub- 
lished work. For example, h258 was the focus of how 
a large di sk can regrow by 2 = after a low z major 
merger in IGovernato et all (j2Q09f ). The results in this 
paper are based on the same initial conditions as those 
previous studies, but simulated at a higher resolution 
that allows for the inclusion of new physics, discussed 
below. The spline force softening of the high resolution 
regions of both h277 and h258 is 174 pc. High resolu- 
tion dark matter particles have masses of 1.3x10^ Af©, 
while gas particles start with 2.7 xlO'^ Mq. Star parti- 
cles are born with 30% of the mass of their parent gas 
particle (i.e., a maximum initial mass of 8100 Mq), and 
lose mass through SNe and stellar winds. Each of these 
galaxies has roughly 5 million DM particles within the 
virial radius at z=0, and more than 14 million particles 
(dark matter 4- gas -f- stars) total. 

The mass and force resolution o f these runs is com - 
parablc to the "Eris" simulation (jGuedes et al.|[20Tll ). 
one of the highest resolution N-body-|-SPH cosmologi- 
cal simulations of a Milky Way-mass galaxy run to date. 
However, unlike Eris, we take advantage of the increased 
simulation fidelity that sufRciently resolves the high den- 
sity regions where stars form [p ~ 100 amu/cm'^) to alter 
the gas cooling and star formation prescription to include 
metal line cooling and H2. The inclusion of metal line 
cooling allows much more gas to cool to the central re- 
gions of the galaxy (Christensen et al., in prep.) and 
become fuel for star formation. However, the star for- 
mation in the new model is tied directly to the local H2 
abundance, which is regulated by the gas metallicity and 
the ability of the ga s to self-shield, in accordance with ob- 
servational results (ILerov et all 12 008: Bi giel et all 120081: 
IBlanc et al.|[2009tlBigiel et al.ll2010;.Schruba et al.ll201lD . 
The H2 prescription imple r nented here is described in 
detail in iChristensen et al.l (|2012[ ). Briefly, it includes 
both gas-phase and dust grain formation of H2 , and de- 
struction of H2 via photodissociation by Lyman- Werner 



radiation from nearby stellar populations. Before the 
implementation of the H2 model, gas particles were re- 
quired to be above a set density threshold and below a 
set temperature before they could form stars. Tying star 
formation to molecular gas eliminates the need to set a 
density threshold above which stars can form, as it en- 
sures that stars form at high densities. A gas particle 
must be colder than lOOOK before it can spawn a star 
particle. 

Other aspects of the cooling physics and energy feed- 
back fro m supernovae rem ain unchanged in these simu- 
lations. IShen et all (|2010D describes the primordial and 
metal-line cooling used, as well as the diffusion of met- 
als that captures the effect of the turbulent interstellar 
medium on metal mixing. A uniform UV background 
turns on at z = 9, mimic king cosmic reionization follow- 
ing a modified version of lHaardt fc Madaul (I200I . Star 
particles are born w ith a Kroupa initial mass function 
(|Kroupa et allll993[) . Each SN deposits 10^^ ergs of ther- 
mal enfirg;5;_withina^ radius calculated follow- 
ing [0itriker2|^cKei (|1988f ). with cooling turned off in 
the affected region for a time that corresponds to the ex- 
pa nsion phase of t he SN remnant, and desc ribed in detail 
in IStinson et al.l (I2006D . As described in IStinson et al.l 
(I2006D. supernova type l a an d II yields are adopted from 
[Thielemann et al.l (|1986[) and lWooslev fc Weaveil (I1995D 
respec tively, and implemented following iRaiteri et al.l 
(jl996[) . This prescription was shown to reproduce the 
observed stellar mass -- metallicity relatio n for galaxies as 
a fun ction of redshift (jBrooks et al.i,2007; iMaiolino et al.l 
[2001 . 

We demonstrate in this paper that including baryons 
in these simulations can dramatically alter the satellite 
evolution in comparison to the dark matter only case. 
Hence, we also have DM-only runs for the two halos de- 
scribed above. The force resolution and z = properties 
(i.e., number of DM particles and halo mass) of these 
DM-only halos are identical to the SPH runs, though 
the mass of the dark matter particles in the SPH runs is 
lower by a factor of (1 — /{,), where is the cosmic baryon 
ratio, ilb/fi,„, and is 0.175 for the adopted cosmology. 

2.1. Satellite Selection & Luminosity Function 

Because the main purpose of this work is to study 
the properties of the satellites of Milky Way (MW)-mass 
galaxies, we explain our method for selecting satellites 
in this section. Halos and their subhalos are identified 
at eac h output time step (every 320 million years) using 
AHF ()Knollmann fc Kne"bi[200l iGiU et al.([200l . For 
field galaxies, AHF adopts overdensities with res pect to 
the cr itical density as a function of redshift from iGrossI 
(jl997[ ). At z=0 the virial radius of h258 and h277 are 
240 kpc and 230 kpc, respectively. 

We first find all of the subhalos that are within the z =0 
virial radii of the two primary halos described above. We 
then reduce that list to only those subhalos that are lu- 
minous, with a minimum of 10 star particles. This list 
of subhalos is likely to include some subhalos that have 
undergone a l arge a mount of dark matter stripping, but 
IBrooks et al.l (|2007f ) showed that the star formation his- 
tories (SFHs), and hence stellar masses, of these simula- 
tions converges only when the halo has more than -^3500 
DM particles. Hence, we then trace back the primary 
progenitor halo for each of the z = subhalos. To do 
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TABLE 1 
Satellite Properties 



Satellite ID 


Mv 


Mhi[Mq] 


Mhi[Mq] 


^infall 






2 = 


z = 










(1) 


(2) 


(3) 


(4) 


h258, 


satl 


— 13.7 


3.5e5 


1.5e7 


0.6 


li258, 


sat2 


— 12 8 








h258, 


sat3 


-12.6 


737.5" 


6.5e6 


2.1 


h258, 


sat4 


1 Q Q 

— io.o 


U. 


o.4eD 


1 Q 
i.O 


h258, 


sat5 


— io.U 


U. 


0. ( e ( 


O 1 


h258, 


sat6 


— il.o 


1 o^c^ 
i.ZCD 


Z.oeo 


1 n 


h258, 


sat7 


1 O 1 


U. 


n 1 -ic; 
y. ico 


1 o 


h258, 


satS 


— y.u 


U. 


U 


1 c; 
I.O 


h258, 


sat9 


— 9.9 


0. 


2.2e5 


1.8 


h258, 


satlO 


—8.8 


0. 


Q 


1.3 


h258, 


sat 11 


-10.7 


0. 


1.4e6 


1.0 


li258, 


sat 12 


-9.5 


0. 


5.4e4 


2.1 


h258, 


satis 


-10.1 


0. 


7.3e3 


1.2 


h277, 


satl 


-14.7 


3.2c7 


2.2e7 


0.02 


li277, 


sat2 


-13.9 


3.8e7 


3.7e7 


0.2 


h277, 


satS 


-14.3 


5.1e6 


3.8e7 


1.5 


h277, 


sat4 


-13.4 


1.3e6 


3.0e6 


0.9 


h277, 


sat5 


-10.4 


0. 





1.4 


h277, 


sat6 


-10.6 


0. 


9.0e6 


1.6 


h277, 


sat7 


-9.9 


0. 


1.4e5 


1.6 




Note. — Column (1) V-band magnitude, calculated based 
on the age and nietallieity of the star particles and adopt- 
ing the Starburst99 stella r population synthes i s mod els of 
iLeitherer et a O l[^9g f) and IVazquez Leithererl l |2005f) for a 
Kroupa IMF ||Kroupa|200lD . Column (2) is the mass in HI gas 
in the satellite at z—0. Column (3) is the mass in HI gas in 
the satellite at infall. The redshift of infall is listed in column 
(4). 

" Given the radial distance of this object from the center of 
tlie main halo that it orbits, this mass of HI would remain 
undetected by current HI observations. We therefore consider 
it a dSph analog. 

this, we successively identify the halo at each higher z 
step that contains the most DM mass of the lower z halo. 
We verified that all of our z = luminous satellites have 
had more than 7000 DM particles within their virial ra- 
dius at some point in the halo's history, ensuring that 
their stellar masses are robust to resolution effects. We 
exclude bright, Magellanic-like satellites with My < —15 
from this analysis, which removes from each simulation 
one bright satellite. However, we make no cut on satellite 
morphology at z = 0, so that the final list of subhalos 
includes both gas-free and gas-rich satellites. 

Table 1 lists some of the properties of the satellites 
studied in this paper. We note that h258 has 13 lumi- 
nous satellites and h277 has 7 satellites, despite their 
similar halo masses. This demonstrates that the number 
of satellites (and particularly the number of highest mass 
satellites) is very stochastic at a fixed pa rent halo mass 
(IVera-Ciro et alJlMl ISawala et al.l[2012h . 

The orbital evolution of every satellite was traced with 
respect to the parent halo. Although some satellites may 
enter the virial radius of the parent halo and exit again, 
we identify the infall redshift as the time at which the 
satellite first enters the parent halo's virial radius. The 
infall times are listed in Table 1, a nd in all cases the 
infall redshift is at z < 3 (see also IZentner &: Bullocl3 
I2003t iGeen et al.ll2012D . We identified the halos in the 
DM-only runs corresponding to each of the dwarf satel- 
lites in the SPH runs by first identifying the halos with 
the best position matches at z=S (prior to infall and any 



Fig. 1. — The satellite luminosity function for our two simulated 
galaxies compared to the Milky Way and M31 satellite luminosity 
functions. The simulated luminosity functions contain all satellites 
listed in Table 1. 



stripping). We then confirmed the matches by ensuring 
that the orbital histories of the matched halos around 
the primary galaxy were comparable. For three of our 
satellites, a good match was not found, and these three 
are excluded from any direct comparison of matched sub- 
halos in the remainder of the paper. 

Figure [1] shows the satellite luminosity functions for 
our two simulated galaxies compared to that of the Milky 
Way and M31. The simulation luminosity functions in- 
clude both gas-free and gas-rich satellites. The MW and 
M31 luminosity functions include all satellites fainter 
than Aly = —15, including Sagittarius a nd Canis Ma- 
j or in the case of the MW, as compiled in iMcConnachid 
()2012|) . Our simulated gas- free satellite sample spans 
the luminosity range of the Milky Way's classical dwarf 
spheroidals, —13.3 < My < —8.8, from Fornax to 
Canes Venatici I (M y = —13.4 and —8.6, respectively, 
lMcConnachi3[20Ta ). 

We stress that no attempt was made to explicitly 
match the classical dSph luminosity range, although our 
final sample does. Instead, the deposition of SN en- 
ergy combined with H2-based star formation was im- 
plemented to reproduce the stellar mass of the parent 
halo at the given halo mass. While past simulations 
have overproduced stellar mass at a given halo mass 
(iZolotov et al.|[200i iGuo et al.llMoHSawala et al.l l2011l: 
iBrooks et al.ll2011l: lLeitneiil2012r it has been suggested 
that restricting star formation to depend on the mass 
in molecular gas (rather than total gas mass) will alle- 
viate this problem, particularly at high z where metal- 
licities are low and formation of H2 on dust grains 
is reduced (iRobertson fc Kraytsovl l2008t iGnedin et al.l 



200llGnedin fc Kravtsovl2010HKrumholz fc Dekelll2011l : 



Feldmann et al.l 120111: iKuhlen et al.l |2012[) . Munshi et 
al. (in preparation) demonstrate that the parent halos 
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stellar mass to 



halo mass relation ( 


Guo et al.ll2010l:lBehroozi et al.ll2010l: 


Moster et al.ll2010l: 


Comov & WcchsW2009l;lYanE et alj 


201lHNeistein et a 


Jl201lHLcautliaud ct al. 2012). Also, 



since star formation becomes metallicity dependent in 
this H2-based scheme, it lowers the star formation effi- 
ciency in the low mass, low metallicity satellites stud- 
ied in this paper, preventing them from overproducing 
stars. iGovernato et al.l ()2012|) previously demonstrated 
that the stellar masses and total masses at < 1 kpc for 
simulated field galaxies in the same stellar mass range as 
the satellite galaxies in this paper are in excellent agree- 
ment with the observational data. 

It is possible that some of the gas-rich satellites in- 
cluded in this work may be unphysically gas-rich, because 
SPH docs not accurately model small-scale hydrodynam- 
ical instabilitie s that can lead to enhanced gas strippin g 
from sate llites (lAgertz et al. l 120071: IMitchell et al.l 120091 ). 
However, IWeisz et al.l ( 2011 ) have recently shown that 
dSphs with Mb < — 10 have extended SFHs, suggesting 
that the majority of observed gas-free satellites at 2 = 
had gas until recently. Thus, despite the fact that some of 
our satellites in the luminosity range of the MW's dSphs 
have gas at 2; 0, their SFHs arc consistent with those 
observed. 

3. REDSHIFT RESULTS 

In this Section, we present the density and circular ve- 
locity profiles of satellites of MW-massed disk galaxies 
at present day. We begin in Figure [U where we com- 
pare the circular velocity, Wc, curves of all the satellites 
in the SPH and DM-only runs. The impact of baryonic 
physics on the circular velocity curves of the simulated 
satellites is immediately clear in this figure. Despite be- 
ing run with the same initial conditions, the subhalos in 
the SPH runs all have lower circular velocities at 1 kpc 
than satellites in the DM-only runs. Since measuring Vc 
at a given radius of a galaxy is equivalent to probing the 
enclosed total mass at that radius, this figure highlights 
that SPH satellites contain less mass within their inner 
1 kpc than their DM-only counterparts. 

In order to study the effect that star formation has had 
on the internal structure of the satellites, we have divided 
the satellite sample by stellar mass at z = 0. "Luminous" 
satellites have stellar mass greater than lO'^Mo, while 
"low-luminosity" satellites arc those with stellar masses 
less than 1O^M0. We further subdivide the luminous 
satellites into categories of gas-rich and gas-free. We do 
this in order to distinguish between simulated satellites 
that would match the MW's dSph population, which are 
all gas-free dwarfs, and simulated satellites that more 
closely resemble the Magellanic clouds, which are gas- 
rich. All of the satellites categorized as "low luminosity" 
are gas-free. 

In Figure [3] we show the dark matter density pro- 
files, p{r), for a subsample of the SPH satellites (black 
solid lines) and their DM-only counterparts (blue dashed 
lines) . Each panel in this figure shows satellites from one 
of the three categories we have defined - the three most 
luminous gas-rich satellites are shown in the left panel, 
the three most luminous and gas-free galaxies in the mid- 
dle panel, and the three least luminous satellites in the 
right panel. Figure |3] shows that the most luminous and 
gas rich satellites in the SPH runs are less dense and have 
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Fig. 2. — The circular velocity, Va, profiles of satellites in DM-only 
simulations (blue dashed lines) and baryon -|- DM (SPH) simula- 
tions (black solid lines), run from the same set of initial conditions. 
The SPH satellites have lower circular velocities at 1 kpc than satel- 
lites in the DM-only runs. 

flatter inner density profiles than their DM-only coun- 
terparts. While the luminous gas-free satellites (middle 
panel) do not appear strongly cored, a direct comparison 
between the SPH and DM-only matches shows that the 
SPH runs do have a more shallow inner profile. Impor- 
tantly, these gas-free, luminous satellites have dramati- 
cally lower densities overall than their DM-only counter- 
parts. At the lowest luminosity end, however, the satel- 
lites in the SPH and DM-only runs tend to have compa- 
rable density profiles (slope and normalization). Hence, 
the process that lowers the density in our SPH satellites 
is more effective at the high stellar mass end than at the 
low stellar mass end. We demonstrate that this process 
is related to the SFH of each satellite in the next sec- 
tion. We note that because the DM particle masses in 
the SPH run are lower than in the DM-only run by the 
cosmic baryon factor, fbar (i-e., these particles have been 
split into DM and gas in the initial conditions), we have 
reduced the DM-only densities by fbar in order to make 
a direct comparison in this figure. 

Figures [2] and [3] demonstrate that, even at moderate 
dSph galaxy luminosities (My < —12, similar to Leo 
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Fig. 3. — The DM density profiles of SPH satellites (solid black lines) and their DM-only counterparts (blue dashed lines) at ^ = 0. The 
left panel shows the three most luminous SPH satellites, which are also still gas-rich at ^ = 0. The middle panel shows the three most 
luminous, gas-free satellites, and the right panel shows the three least luminous SPH satellites. Luminous SPH satellites have significantly 
shallower central density profiles than DM-only satellites, while low-luminosity SPH satellites retain central density cusps similar to their 
DM-only counterparts. 



I or And II), baryonic processes result in an expanded 
and shallower central dark matter distribution than pre- 
dicted by DM-only simulations. Even when the corecfl 
density profiles of our simulated dwarfs do not have a 
flat slope (7 = 0), the absolute values of density in the 
central regions are still dramatically reduced compared 
to expectations from a DM-only simulation. We con- 
clude that DM-only simulations, or models based on an 
assumption of adiabatic contraction, make physically in- 
correct predictions for the central (< 1 kpc) masses of 
dwarf galaxies more luminous than Leo I {My = —11.9). 
Such models should therefore be used with caution when 
used to interprete the observed dynamics of dwarf satel- 
lites in a cosmological context. 

The general conclusions disc ussed in this pa,per ar e in- 
dependent of resolution effects. iSpringel et al.l (|2008[ ) ex- 
amined the convergence of density profiles in very high 
resolution DM-only simulations (the Aquarius simula- 
tions), and found that the density converged at 4 ~ 6x 
the force softening, e. We verified, using a lower resolu- 
tion DM-only run with 1.5 x larger force softening, that 
the density slopes in these simulations converge by 4e, or 
700pc. Using the same lower resolution simulation, we 
found that Vc is 90% converged by 1 kpc. The Vc con- 
vergence radius is larger than for density, because Vc is a 
cumulative quantity. 

In the remainder of the paper we will often compare 
Vc values in the SPH and DM runs at 1 kpc. As will be 
seen in the next section, the Vc values at 1 kpc show the 
dramatic effect that baryonic physics has on the SPH 
runs, but avoids biasing this value due to convergence 
issues at smaller radii. More importantly, the Vc val- 
ues at 1 kpc between the SPH and DM-only runs are in 
excellent agreement in our low-mass subhalos for which 
baryons do not dramatically alter the evolution. The 
convergence of SPH and DM-only results in these lower 
mass halos (which, by definition, contain fewer particles 
and are less resolved than their high mass counterparts) 
demonstrates that there are no spurious numerical ef- 
fects introduced by the lower mass baryonic particles in 

® In the rest of this paper, we will use the term "cored" to refer to 
all slopes shallower than predicted by DM-only simulations alone. 



the SPH runs. 

4. MASS DEPENDENT EVOLUTION OF SATELLITES 
WITH BARYONS 

In this Section, we study the evolution of the satellites 
to understand the processes that lead to the lower con- 
centration of mass at z =0 in the SPH satellites. We 
first focus on the evolution at high redshift. and demon- 
strate that DM core creation occurs in the most luminous 
satellites prior to their infall. After infall, we show that 
tidal stripping effects exacerbate the mass discrepancy 
between SPH satellites and their DM-only counterparts. 

4.1. The Impact of Baryons Before Infall 

We now examine the evolution of Vc from high z to in- 
fall for all of the satellites in our sample. The top panel 
of Figure [5] shows the change in the DM contribution to 
Vc at 1 kpc between SPH satellites and their DM-only 
counterparts at infall, as a function of the stellar mass of 
the SPH satellites at infall. It can be seen in this panel 
that satellites that have formed more stars prior to infall 
(those with M^ > lO^M© at infall) under go a significant 
decrease in DM mass interior to 1 kpc, in comparison 
to their matched counterparts in the DM-only runs. The 
DM rotation curves of these luminous satellites arc there- 
fore 2 — 16 km/s lower than those of DM-only satellites. 
(Masses in the DM-only run have been reduced by f^ar 
for a direct comparison.) We conclude that baryonic ef- 
fects lower the central DM densities, and hence lower the 
central DM circular velocity, of massive satellites prior 
to infall. 

The bottom panel of Figure 2] shows the total change 
in Vc at 1 kpc between SPH satellites and their DM- 
only counterparts at infall. Here we find essentially no 
trend in the change of the total Vc of SPH satellites 
with stellar mass. This is because at the luminous end 
(M^ > IO^M^q), baryons contribute significantly to the 
total mass within 1 kpc. However, once the mass con- 
tribution from baryons is subtracted, it is seen that the 
remaining mass in DM is significantly reduced by bary- 
onic effects. We note that the most luminous satellite at 
infall {Mi, ^ 10® Mq) has undergone adiabatic contrac- 
tion, resulting in a lower total Wc(l kpc) than its DM-only 
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Fig. 4. — The difference between Vc at 1 kpc in the SPH and 
DM-only runs at infall, as a function of the stellar mass in the SPH 
satellite at infall. Top panel: The difference in the DM contribution 
to Vc at 1 kpc for matched SPH and DM-only subhalos. Bottom 
panel: The difference in total Vc at Ikpc. 



match^ 

Go vernato et alj (j2012D also found that DM core cre- 
ation varied as a function of mass, for isolated field galax- 
ies. These authors have shown that the creation of DM 
density cores due to SNe driven outflows is common in 
galaxies with M^, > 10^ at z = 0. For lower luminosity 
galaxies, their work finds that SF and its associated feed- 
back arc not efficient enough to have flattened a galaxy's 
steep DM density profile. We find that this mass thresh- 
old above which feedback becomes effective is the same 
in satellite galaxies as well, despite the fact that core 
formation docs not continue after infall. 

Several theoretical models have shown how outflows 
and galactic fountains can lead to the flattening of 
DM density cores in dwarf galaxies fNavarro et all 



1996HRead & Gilmord 


20051: IPontzen & Governatdll20Tl 


Governato et all 12012 


) . Rapid and frequent star 



formation episodes break the adiabatic approxima- 
tion in the central kiloparsec of galaxies, transfer- 
ring energy to the collisionless particles, and result- 
ing in shallow DM density cores (IPontzen fc Governatol 
120121: IQgiva fc Monll2012l ). An essential component to 
such DM-flattening scenarios is the ability to resolve 
the high density clumps where star formation takes 
place, i.e., overdensiti es compa rable to giant molec- 
ular clouds ((Robertson &: Kravtsov. ,20081 : iSaitoh et al.l 



2008t iTasker fc BrvanI [2008t ICeverino fc Klvpinl [200l 
Christensen et al.l 120101: iColm et al.l I2010D . When star 
formation is limited to these high density peaks, energy 
deposited from SNe creates overpressurized regions of hot 
gas, driving outflows of SNe heated gas from the simu- 
lated galaxies. The H2 model adopted in these simula- 
tions includes self-shielding of cold gas, preventing heat- 
ing from photoionization in dense regions with no young 
stars nearby. Self-shielding, combined with a low tem- 
perature cooling model, allows the already cool gas in 
the galaxy to cool even further, making the star forma- 
tion eve n clumpier and the feedback therefo re even more 
efflcient ()Christensen et al.ll2012l:lSusall2008l Christensen 
et al., in prep.). 

When star formation is limited to high density peaks, 
as in this paper, the SNe feedback after a star forming 
event creates regions of overpressurized, hot gas that will 
shut down additional star formation for a period of time 
until the gas can cool and continue to form more stars. 
This process leads to a bursty SFH. Examples of the re- 
sulting star formation histories for galaxies in this paper 
are shown in Figure |SJ Dwarf galaxies in our highest 
luminosity range are massive enough to retain gas for 
an extended period of time, allowing them to undergo 
multiple star forming and gas loss events. Each event 
pushes successively more DM to larger orbits, gradually 
transforming a cuspy DM density proflle into a flatter 
cored proflle. Hence, extended SFHs that allow for mul- 
tiple bursts of SF lead to effective DM core creation. At 
the lower end of the luminosity range of our satellites, 
the halos are less massive, with shallower potential wells. 
A few are low enough in mass that they lose their gas 
relatively early, though most maintain some gas at least 
until infall. The overall gas mass in these low mass sub- 
halos is substantially reduced at early times, due to a 
combina tion of reionization and initial SF that ejects gas 
(see also lSawala et al.llMol) . The remaining gas in the 
shallow potential well then has a difficult time reaching 
the densities required for star formation. The inability 
to continue forming stars in the lower mass, less lumi- 
nous satellites prevents them from having bursts of star 
formation over extended periods, and they thus retain a 
steeper DM density profile. 

If the mass trends described above are accurate, we 
should see these trends reflected concurrently in the 
SFHs of the galaxies and their reduction in central dark 
matter mass. To examine the evolution of the satellites 
as a function of mass, we focus on the progenitors of 
three satellites as examples, drawn from the three differ- 
ent categories represented in Figure O (gas-rich, gas-free 
and most luminous, gas- free and least luminous). Fig- 
ure |5| shows the SFHs of these three representative satel- 
lites. The SFHs, particularly in the luminous satellites, 
are episodic with bursts, rather than smooth or constant. 
In Figure [SI we show the evolution of the DM contribu- 
tion to Vc at 1 kpc in these same three galaxies. Each Vc 
evolution curve is normalized to the Vc value of the DM 
at infall, where the infall times are marked by vertical 
lines, color coded to each galaxy. 

For the luminous satellites (shown in black and red 
in Figures |S] and jB]) , reductions in the DM contribution 
to Vc at 1 kpc can be seen following strong bursts of 
SF in these galaxies. For example, a DM core begins 
to form in the progenitor of the most luminous satellite 
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Fig. 5. — The star formation histories of three SPH sateUites. Note the different range for the axes. At infall, the sateUites in black (left 
panel), red (middle panel), and blue (right panel) have virial masses of 7.9 X 10'^, 2.7 X 10^, 0.6 x lO'' Mq, respectively, and infall times 
of ^ = 0.02, 1.3, 1.5. In the first two panels, the vertical dashed line marks z =infall onto the parent halo. In the final panel, z =infall lies 
beyond the range of the x-axis, which is truncated before infall for clarity. However, the vertical dashed line in the right panel indicates 
the onset of reionization (i.e., when the uniform UV background turns on). 

ization and z = 3. This is evidenced in Figure |6l which 
shows that the DM contribution to Vc in the central 1 kpc 
stays constant through the entire history of this galaxy, 
both prior to and after infall. 

Finally, wc verified that the SFHs of satellites are either 
completely truncated at infall, or SF continues at a low 
rate with no strong bursts of SF that are associated with 
core formation. Thus, we conclude that the dark matter 
core formation occurs prior to infall. The SFHs of the 
galaxies before infall are intimately tied to their DM den- 
sity profiles, with multiple SF bursts re quired to create a 
DM core in the progenitor of a sat ellite (jRead fc Gilmord 
l2005t iPontzen fc GovernatollMl . 

4.2. Evolution After Infall 

In this Section, we focus on the evolution of the cen- 
tral mass of satellites after infall in order to explore 
whether the presence of baryons in galaxies leads to 
diverging tidal evolution between the SPH and DM- 
only cases. In tracing the orbits of our satellite 
sample, we found that some satellites undergo signifi- 
cant mass loss after accretion, with most of the mass 
loss occurring at pericenter passage where tidal effects 
are strongest. This is true in both the SPH and 
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Fig. 6. — The DM contribution to Vc at 1 kpc in the three SPH 
satellites shown in Figure [5] as a function of time. The Vc value at 
all times is normalized to the value at infall. A luminous gaseous 
(at z =0) satellite is in black (solid line), a luminous gas-free (at 
z=0) satellite is in red (dashed lino), and a low luminosity gas- free 
(at z=0) satellite in blue (dotted line). Their respective infall times 
are shown by the vertical lines, with colors and linestyles according 
to galaxy as just described. In the two most luminous satellites, 
decreases in Vc correspond to bursts of SF seen in Figure [5] 

(shown in black) near z = 1, consistent with a large 
number of repeated bursts of SF that begin at this time. 
Figure ini indeed shows a significant decrease in the central 
DM mass associated with core formation beginning at 
z r^ 1. During periods of low SF, baryons have little 
to no impact on the DM structure of these satellites. 
For example, the satellite shown in red does not undergo 
any changes in the contribution of DM to its central Vc 
between z = 3 and infall, which is tied to the weak SF 
of this galaxy during this period. 

The low luminosity, gas-free satellite (shown in blue in 
Figures [5] and [6|) forms most of its stars prior to reioniza- 
tion. The gas supply of this satellite is strongly affected 
by reionization, with SF declining in the ~ 3 Gyr after- 
ward, until gas is completely unbound from this halo at 
z ^ 2. The single, small burst of SF in this galaxy prior 
to reionization was not enough to alter the DM density 
profile, nor was the steady decline in SF between reion- 



other studies (e.g., 


Maver et al.ll2001l 


Dekel et al.ll2003al: 


Havashi et al.ll2003 


jKazantzidis et al. 


2004: Maver et al.l 


20061 Read et all 120061: Bovlan-Kolchin & Ma 2007; 



2010: Kazantzidis et al.. 120111 ). However, most work to 
date on mass loss in satellites has excluded the effects of 
baryons. We demonstrate below that when a satellite's 
orbit is such that it is likely to undergo significant tidal 
stripping, SPH satellites will lose a substantially larger 
fraction of their central mass after infall compared to 
their DM-only counterparts. 

In Figure [7] we show the change in circular velocity, 
between infall and z = 0, at 1 kpc for all simulated satel- 
lites, as a function of the fraction of total mass lost during 
this time. Two important trends emerge in this figure. 
First, SPH satellites tend to undergo a larger decrease 
in their Vc{l kpc) than their DM-only counterparts af- 
ter infall. We find that the relative reduction of Vc be- 
tween SPH and DM-only satellites is highly dependent 
on the orbital parameters of the satellites. For exam- 
ple, satellites accreted more than 6 Gyr ago [zinfaii > 1) 
experience a reduction in their central circular velocities 
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Fig. 7. — The change in Vc at 1 kpc between z =infall and z =0 
for the simulated satellites, as a function of the total mass lost 
during that time. The SPH satellites (in circles) are divided by 
their infall time, and the matched DM-only counterparts are shown 
as blue diamonds. 

that is 11 — 62% more than the reduction in Vc{l kpc) 
experienced by their DM-only counterparts. Likewise, 
satellites on more radial orbits {rp^ri ■ rapo > 8) un- 
dergo a reduction in Vc that is > 35% greater than their 
DM-only matches. On the other hand, SPH satellites 
accreted more recently {zinfaii < 1) or on circular orbits 



{fpe 



' apo 



< 2) have their Vc reduced by only — 13% 



more than their DM-only counterparts. A larger decrease 
in Uc(lkpc) will occur in some of the SPH satellites sim- 
ply from gas loss (e.g., in ram pressure stripping). A 
comparison of the HI masses at infall to the z = HI 
masses in Table 1 demonstrates that significant amounts 
of gas are lost. In the vast majority of the satellites, 
however, the decrease in Vc in the SPH satellite is signif- 
icantly higher than can be accounted for by just loss of 
gas. The second trend seen in Figure [7] is that only SPH 
satellites lose more than than 90% of their total massFl 
Overall, we find that SPH satellites lost 3 — 34% more of 
their virial mass than their DM-only counterparts after 

i nfall. 

iPefiarrubia et al] (120101 hereafter PIG) used N-body 
simulations to examine the role of different DM density 
profiles, as well as the role of a baryonic disk in the host, 
on the tidal evolution of satellites. Their models show 
that presence of a baryonic disk in the primary galaxy 
results in a higher mass-loss rate for satellites at each 
pericentric passage, in comparison to a host with no disk. 
The infiucncc of a baryonic disk is especially strong for 
satellites with shallow DM density profiles, with slopes of 
0.0 < 7 < 0.5, resulting in a higher fraction of mass loss 
at each pericenter passage for the cored satellites. The 
combined effects of the host disk and a satellite's shallow 
DM density profiles in PIO becomes increasingly signifi- 



Subhalos must have lost ~90% 
stars, which are mor e tightly bound, 



of their mass before 

^ _ _ begin to be stripped 

llPeiiarrubia et al.l2008l) . Only three satellites in our sample change 
their stellar mass by more than 20% between infall and z =0, and 
they are the same three that have lost more than 90% of their total 
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Fig. 8. — The Vmax values of the SPH subhalos at infall as a 
function of their stellar mass at infall. The solid line shows the 
relation M* <x vf^^^, normalized at M* = 10^ Mq, indicating how 
stellar mass increases with Vmax in these simulations. 

cant with increasing time after infall and for satellites on 
more eccentric orbits. 

The greater reduction in central circular velocities and 
virial masses we find in our SPH satellites is due in part 
to the presence of a disk in the SPH host galaxies (the 
DM-only runs do not have a disk) , which results in more 
efficient mass loss for the SPH satellites. As noted in PIO, 
the effect of the disk should be even more pronounced in 
the tidal evolution of cored satellites, since the disk can 
dominate the tidal field in the center of cored satellites, 
but not in the tightly bound cuspy halos. The small 
sample size of our satellites and the strong dependence 
of tidal stripping on orbital parameters makes it difficult 
to test this directly in this paper. However, our sample 
does contain two SPH satellites - one with a cored DM 
density profile and one with a cuspy DM density profile 
- that have similar infall times and very similar orbital 
parameters. The Vc at 1 kpc of the cored SPH satellite 
was reduced by 78% after infall, while the cuspy SPH 
satellite underwent only a 44% reduction in its central 

Vc- 

Wc conclude that if a satellite is likely to have under- 
gone significant tidal stripping due to its orbital history, 
the existence of a baryonic disk in the host and a shallow 
inner DM density profile in the satellite will exacerbate 
the amount of mass lost. 

5. AN UPDATE TO THE THEORETICAL MODEL 

When using DM-only theoretical results to interpret 
the MW's dwarf satellite population, it is common to as- 
sociate the DM-only subhalos that were the most massive 
at the time of their formation or accretion with the more 
luminous dSphs (jS ullock et al. 



,_ 2000; Kravtso v et al l 

20041: iGnedin fc Kravtsovl 12 006: Strigari et al. 2007al: 
Koposov et al.ll2009l: iBoviU fc Ricotti,2011 - Simlia et al.l 
2012D . i.e., an abundance matching technique. This as- 



signment is reasonable if we assume that the most mas- 
sive halos will have the highest star formation rates, and 
therefore be the most luminous at z = 0. Figure |8] shows 
that our simulated halos follow a trend of increasing stel- 
lar mass with halo mass (represented by Vmax) prior to 
infall, so that the most luminous galaxies at z = in- 
deed correspond to the most massive galaxies at infall. 
We verified that the matched DM-only counterparts can 
be assigned to the most luminous satellites at infall in 
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the SPH run. In the se simulations, M, oc M^^^ (see 
iGovernato et al.|[2012[) . Because v^ax roughly scales as 

(e.g.. IKlvT)in"eFaII[20lTI) . cx vf^^^. This rela- 
tion is shown as the solid line in Figure [HJ normalized 
at M* = 10^ Mq. We note, however, that tidal strip- 
ping can substantially reduce the mass of some halos by 
z = 0, introducing scatter into the tight relation seen in 
Figure [S] at infall. 

While the assumption that stellar mass increases with 
halo mass is valid for these simulations, it cannot be 
assumed that the central dark matter mass distribution 
of the DM-only subhalos is the same as that of luminous 
subhalos. Despite the fact that the virial masses at infall 
arc similar in the SPH and DM-only runs, the inner DM 
orbits have been expanded in the SPH runs (but, prior 
to infall, arc not removed from the halo). As this paper 
has shown, the DM mass in the inner regions is lowered 
for those satellites with M* > IG^Af© due to the effects 
of feedback prior to infall. After infall, SPH satellites are 
also prone to significantly more tidal stripping than their 
DM-only counterparts. Figure |9] summarizes our results 
on the combined impact of DM core-creation and tidal 
stripping on the internal dynamics of satellites. This 
figure shows the difference in Vc at 1 kpc at 2; = between 
SPH and DM-only satellites, as a function of Vmax at 
infall of the DM-only satellites. The top panel of this 
Figure shows the change in the DM contribution to Vc 
(again, the DM-only masses have been reduced by fbar 
for a direct comparison), while the bottom panel shows 
the change in the total Vc at 1 kpc. We note that we 
exclude from the plot one data point that has a change 
in total Vc 25 km/s. This dramatic change is due to the 
fact that the SPH satellite has recently passed directly 
through the disk of its parent galaxy, stripping it more 
substantially than its DM-only counterpart (since there 
is no disk in the DM-only run). 

A well defined trend for the change in DM mass interior 
to 1 kpc between SPH and DM-only runs can be seen in 
the top panel of Figure |9l In order to quantify the change 
in DM mass that baryonic processes combined with tidal 
stripping induce in the central regions of satellites com- 
pared to the DM-only case, we fit a linear regression to 
the data points in the range 20 < Vmax < 50 km/s. The 
resulting fit, shown in the top panel of Figure 19] as the 
dashed line, is A(wc, Ikpc) = 0.2vm.ax,DM-oniy - 0.26. 

The bottom panel of Figure [9] shows that the change 
in total Vc at 1 kpc between SPH and DM-only runs is 
slightly lower for gaseous satellites (open circles) than 
those satellites that are gas-free. Like the bottom panel 
of Figure m this is because baryons contribute signifi- 
cantly to the mass in the interior Ikpc. While there is 
not much trend at infall visible in the bottom panel of 
Figure m a trend is beginning to emerge at z = in the 
bottom panel of Figure |9l That is, as gas is removed 
from these systems after infall, the underlying difference 
in the DM central masses in the SPH and DM-only runs 
emerges. Note that it is likely that the gas-rich satel- 
lites in this work are artificially gas rich due to ineffi- 
cient stripping of gas in SPH simulations. The trend in 
the bottom panel of Figure El therefore, would be even 
stronger if SPH did not suffer from this numerical effect. 

Wc advocate that DM-only results apply the above lin- 
ear relation to those halos that are accreted onto Milky 
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Fig. 9. — The difference in Vc at 1 kpc at z = between the SPH 
and DM-only counterparts, as a function of Vmax of the DM-only 
satellite at infall. Top panel: The difference in the DM contribution 
to Vc at 1 kpc for matched SPH and DM-only subhalos. Bottom 
panel: The difference in total Vc at Ikpc. 

Way-mass halos with Vmax > 20 km/s. Failure to ac- 
count for these effects will over predict the central mass 
of luminous satellites such as those around the Milky 
Way and M31. This is likely the source of th e tension dis- 
cussed in lBovlan-Kolchin et al.l (j201lL[20T^ between the 
densities in the MW's classical dSph population and the 
DM-only Aquarius subhalos. Those studies emphasized 
that the observed MW dSphs appear to be less dense 
than those found in a CDM DM-only simulation. In- 
cluding the affects of baryonic physics reduces, and may 
even completely alleviate, this discrepancy. We explore 
the observational consequences of our updated model fur- 
ther in a companion paper (Brooks & Zolotov 2012). 

6. CONCLUSION & DISCUSSION 

We have demonstrated that the inclusion of baryonic 
physics can dramatically alter the evolution of bright 
satellite galaxies around Milky Way-massed galaxies. 
The population of satellites studied here have F-band 
magnitudes consistent with the range of the Milky Way's 
classical dSphs, —15 < My < —8. Our sample contains 
both gas-frec dSph analogs and massive gas-rich dwarfs. 
By directly comparing the internal properties of satellites 
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simulated with gas hydrodynamics (SPH) to the same 
satellites in DM-only simulations, we have demonstrated 
the impact of baryons across a range of masses on satel- 
lite evolution. Our main results are summarized below. 

• Before infall, the progenitors of luminous satellites 
{My,r > 10^ Mq,M^ > IO^Mq) undergo rapid and 
frequent bursts of star formation. The associated 
SNe feedback from these SF episodes results in 
structural changes to the central mass distribution 
of these dwarf galaxies, resulting in reduced DM 
densities and shallower inner DM density profiles 
than DM-only galaxies. 

• The progenitors of lower luminosity satellites 
(M^ < lO^Af© at infall) reduce their gas con- 
tent, and hence become less efficient at forming 
stars, earlier than more massive dwarfs. This 
gas loss is partly due to heating by the uniform 
UV background, as well as subsequent gas loss 
in early star forming/feedback events, preventing 
them from having multiple strong bursts of star 
formation that lead to DM core creation. Low lu- 
minosity satellites, therefore, tend to retain steep 
DM density profiles that are comparable to DM- 
only runs. 

• Once accreted, SPH satellites experience more 
mass loss due to tidal stripping than DM-only satel- 
lites, the amount of which is dependent on infall 
time and orbit. While both SPH and DM-only 
satellites are affected by tidal stripping, the pres- 
ence of a baryonic disk in the SPH runs results in a 
greater reduction in the central Vc in SPH satellites. 
The influence of a baryonic disk is especially strong 
for satellites with shallow DM density profiles. We 
find that SPH satellites with Zinfaii > 1 experience 
a reduction in their central circular velocities that 
is 11 — 62% more than the reduction in Vc at 1 kpc 
experienced by their DM-only counterparts. 

• By z = 0, the combined effects of DM core creation 
and enhanced tidal stripping for luminous satel- 
lites results in a significant discrepancy between 
the circular velocity profiles of SPH and DM-only 
satellites. We find that the Vc at 1 kpc predicted 
for satellites by DM-only simulations should be re- 
duced by A{Vc,lkpc) 0.2Vmax,DM-only - 0.26, 

for satellites with 20 < v„iax < 50 km/s at infall. 

High resolution that allows simulators to limit star for- 
mation to high density peaks is an essential requirement 
to re produce the bary o nic eff ects in this paper. How - 
ever. iGovernato et aLl (|2010( ) and iGuedes et all (|2011l ) 
showed that, even at high resolution, if star formation 
is allowed to occur diffusely across the disk, no large 
scale outflows are generated. Restricting star formation 
to high density peaks instead leads to overpressurized 
regions of hot gas when stars go SNe, leading to out- 
flows. These overpressurized regions expand faster than 
the local dynamical time. When this occurs in the cen- 
tral ^1 kpc, the potential flattens as this hot gas ex- 
pands, lBadmg_to_a£_jTrevCT of the DM 
orbits (jPontzen fc Governatci|[20T^ . Restricting star for- 
mation to high density peaks is comparable to allowing 



stars to only form in giant molecular clouds, rather than 
across the entire disk at any given time. The simulations 
used in this work are the flrst at these high reso lutions 
to include metal line cooling (jShen et al.l l201Ci[ ) and a 
presciption for self-shielding of cold gas, allowing star 
formation to be tied to the shiel ded regions where H2 
can form (jChristensen et a"Lll2012|) . 

It is important to note that simulations that do not 
resolve the effect of feedback at high densities will be un- 
able to reproduce the results of this paper. However, 
the feedback model employed in this paper has been 
shown to match the observed mass - me tallicity rela- 
tion for galaxies as a func tion of redshift (jBrooks et al.l 
l2007tlMaiolino et al.|[2008l ). the baryonic TuUy-Fisher re- 
lationship (Christensen et al. in prep.), the size - lu - 
minosity relation of galaxy disks (jBrooks et al.l 120 lit ) , 
the z = stellar mass to halo mass relation (Munshi 
et al., in prep), and the central mass as a function of 
stellar mass for galaxies in t he luminosity range in this 
paper (jGovernato et al.ll2012D . This large number of suc- 
cesses in matching the observed scaling relations of galax- 
ies lends credence to the particular feedback model em- 
ployed in this paper to study satellite galaxies. 

In addition to reproducing the above scaling relations, 
the star formation and feedback model used in this work 
has been used to simu late bulgeless dwarf disk galaxies 
(|Governato et al.||2010D . Importantly, the processes that 
lead to bulgeless galaxies also transform cuspy DM den- 
sity profiles into cored profiles, leading these simulations 
to match the central dark matter densities d erived by 
the THINGS and Little THINGS surveys (|0h et all 
I2OIII : IGovernato et al.ll2012l ). In other words, the simula- 
tions used in this work have been shown to reconcile the 
cusp/core problem in CDM. In this paper, we have shown 
that this same model can alleviate the tension between 
the dense, massive satellites predicted by CDM with 
the o bservations of lower density , luminous dSph satel- 
lites (IBovlan-Kolchin et all lMa IWoff fc Bullock 12011 
IHavashi fc Chiball2012n . We note that the processes de- 
scribed in this paper also act to reduce the overall num- 
ber of massive subhalos that exist at z = 0, potentially 
solving the missing satellites problem in CDM. Hence, 
it remains possible to resolve the small scale problems 
of CDM with a proper model for baryonic physics, and 
without invoking exotic forms of DM. 

The results presented here show that using DM-only 
CDM simulations to study the internal dynamics of lu- 
minous satellites will lead to erroneous results. While it 
is safe to assign the most luminous satellites to the orig- 
inally most massive halos, those massive halos will expe- 
rience evolution that CDM DM-only runs don't account 
for. In Brooks & Zolotov (2012). we address how the 
model presented here affects the interpretation of kine- 
matic observations of the Milky Way's dSph population. 
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